Complexity and Fragility in Ecological Networks 
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A detailed analysis of three species-rich ecosystem food webs has shown that they display scale- 
free distributions of connections. Such graphs of interaction are in fact shared by a number of 
biological and technological networks, which have been shown to display a very high homeostasis 
against random removals of nodes. Here we analyse the response of these ecological graphs to both 
random and selective perturbations (directed to most connected species). Our results suggest that 
ecological networks are extremely robust against random removal but very fragile when selective 
attacks are used. These observations can have important consequences for biodiversity dynamics 
and conservation issues, current estimations of extinction rates and the relevance and definition of 
keystone species. 



I. INTRODUCTION 



Ecological research has widely demonstrated that com- 
munity fragility is far from being understood. Issues as 
which species might be considered as specially relevant 
because of their strong effects on the community have 
lead to a heated debate since Paine's definition of key- 
stone species (Jordan et al.,1999). Despite this and other 
discussed topics, it is commonly accepted that commu- 
nity fragility is related to how ecological communities 
are structured, specifically to how trophic links are dis- 
tributed throughout the community (May 1974; Pimm 
1991). But both the scarcity of high-quality data (Po- 
lis, 1991; Cohen et al., 1993; Williamson and Martinez, 
) and the lack of methods suitable for a detailed anal- 
ysis of the complexity of food web organization (Cohen 
et al., 1993) leads to a lack of an unified picture of com- 
munity fragility. A number of questions emerge from 
these studies: How are dynamic and static (graph-level) 
properties related?; How dependent is ecosystem fragility 
from graph architecture? 

Recently, there has been an increasing interest in the 
organization of complex networks. These networks go 
from technological ones (Watts and Strogatz, 1998; Al- 
bert et al, 1999), to neural (Watts and Strogatz 1998; 
Amaral et al., 2000) or metabolic networks (Jeong et al 
2000; Wagner and Fell, 2000). All these networks can be 
represented as a graph consisting of a set of nodes and 
the links connecting them. 

Such complex networks share some topological fea- 
tures, as the so-called "small world" (SW) behavior 
(Watts and Strogatz 1998; Newman, 2000). Some of 
these webs also exhibit scale-free (SF) distributions of 
links. Specifically, the frequency of nodes with k con- 
nections follows a power law distribution P(k) w fc -7 , 
where most units are connected with few nodes and very 
few nodes are highly connected. 




FIG. 1. Schematic representation of an ecological graph: 
few species have many links and many have just one or 
two. This leads to fat-tailed (scale free) distributions of 
connections. This particular example follows a power law 
P{k) fa fc -1,2 . This picture reflects the most common fea- 
tures observed in webs such as those reported here. 

Networks exhibiting SW properties and SF distribu- 
tions of connections present a characteristic response to 
the successive removal of their nodes, related to the way 
removals occur (Albert et al 2000; Jeong et al 2000). 
When nodes are removed at random, the network ex- 
hibit high homeostasis. By contrast, if most-connected 
nodes are successively eliminated, the structure of the 
network reveals an intrinsic fragility that eventually leads 
to a breaking into many small subgraphs. This behaviour 
is not shared by other networks, such as purely random 
ones, where P(k) is Poissonian. It has been demonstrated 
that random networks are equally fragile to the way re- 
moval of nodes is produced. 

The surprising and general nature of these results im- 
mediately suggests their application to ecological net- 
works (figure 1), which have been recently shown to dis- 
play SF behavior (Montoya and Sole, 2000). Here we 
examine the possible consequences for ecosystem stabil- 
ity against different types of species loss. As we will 
see, these networks display the robustness expected for 
long-tailed distributions of connections but also a high 
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fragility against selective species removal. 



II. FOOD WEBS ANALYSED 

Because of the limitations of the available data in terms 
of both taxonomic resolution and size (Polis, 1991; Cohen 
et al., 1993; Williams and Martinez, 2000) our study is 
limited to the three richest and best-described food webs 
available in the ecological literature (figure 2). These are: 
Ythan estuary web (Huxman et al. 1996), Silwood Park 
web (Memmott et al. 2000) and Little Rock Lake web 
(Martinez 1991). 
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parasitoids, 5 omnivores, 19 herbivores and one plant. 

Finally, Little Rock Lake food web corresponds to a 
small lake. It is the largest of the three webs analysed 
here (N s = 182), although it has less taxonomic resolu- 
tion. Only 31% of nodes are real species, being most of 
the nodes genera-level (63%) and the rest corresponding 
to higher taxa. Here < k >= 26.05. 

By using these webs, we have a diverse representation 
of habitats: one food web from a terrestrial habitat (Sil- 
wood park), a freshwater habitat (Little Rock lake) and 
an interface environment (Ythan estuary). All of them 
have some features in common. In a previous study, 
we have shown that these networks display small-world 
properties (Montoya and Sole, 2000). In a graph with 
a small world (SW) topology, nodes are highly clustered 
yet the path length between them is small. In this sense a 
SW stands for a network whose topology is placed some- 
where between a regular and a totally random distribu- 
tion of connections. These networks display a number 
of surprising features and have been suggested to be of 
great relevance in different biological contexts (Watts and 
Strogatz, 1998; Jeong et al., 2000; Lago-Fernandez et al, 
2000). 

In general, SW nets have been shown to provide fast re- 
sponses to perturbations and thus provide a great source 
of homeostasis. However, two of these networks also dis- 
play the power law distribution P{k) w fc~ 7 (figure 2) 
thus belonging to the SF class of networks. Little Rock 
web is also fat-tailed but it displays deviations from the 
power law due to low taxonomic resolution (Montoya and 
Sole, 2000). For food webs displaying fat-tailed distribu- 
tions, perturbations can have unexpected consequences, 
which are explored in the next section. 



FIG. 2. Histograms showing the distribution of links (k) 
per number of species (P(k)) for the three networks analysed 
here (black bars): (A) Ythan estuary, (B) Silwood web and 
(C) Little Rock lake. Webs A and B are shown together with 
the best power-law fit (white bars, see table I for details). 
Little Rock web, shown in (C) displays several bumps (two of 
them indicated) due to low-level taxonomic resolution (Mon- 
toya and Sole, 2000). 

Ythan estuary food web has N s — 134 species, be- 
ing the second largest documented web in the UK. Most 
nodes correspond to real species (88%) while the rest in- 
volve lower taxonomic resolution (all the species of Aca- 
rina or of Brown Algae are lumped together in the same 
node). It is one of the most reviewed food webs through 
ecological literature and the average number of links per 
species is < k >= 8.7. 

Silwood park food web is a very detailed sub-web (all 
nodes but one are real species) of those species associated 
with the Scotch broom Cytisus scoparius in a field site of 
97 hectares. The average number of links per species is 
< k >— 4.75. It includes 154 species: 60 predators, 66 



III. RESPONSE TO SPECIES REMOVAL 

We have simulated two kinds of species removals: ran- 
dom and directed. These correspond to removal of an ar- 
bitrary or the highest connected node, respectively Pre- 
vious studies have shown that the eventual effect of re- 
moval is network fragmentation, which takes place in very 
different ways depending on the type of removal used 
(Albert et al., 2000; Jeong et al., 2000). Community 
fragility has been measured in different ways in relation 
to the fraction of species / that have been already re- 
moved (Figure 3). We have measured the fraction of 
species contained in the largest species cluster S for each 
/; the average size of the rest of the species clusters < s > 
as the food web is being fragmented; and the fraction of 
species becoming isolated due to removals of other species 
on whom their survival depends, which are known in the 
literature as secondary extinctions, and can be used as 
measure of extinction rate (Pimm, 1991). Here f c indi- 
cates the fraction of removed species at which the web 
becomes fragmented into many small sub-webs (figures 
3a-c). 
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FIG. 3. Response of food webs under random (circles) and 
directed (squares) species removals, (a-c): relative size of the 
largest species cluster S (open symbols) and average size of 
the rest of the species clusters < s > (filled symbols) in rela- 
tion to the fraction / of removed species for Ythan estuary, 
Silwood Park web and Little Rock lake food web, respectively. 
Critical thresholds f c are indicated. In (a) and (b), continu- 
ous line consider a species cluster that appear quickly. This 
does not rest robustness to the observed trend, because it is 
similar to what happens for < s > and S near f c if that clus- 
ter is not considered in the calculations, (d-f): extinction rate 
(fraction of secondary extinctions) as a function of / for the 
webs ordered as before. 
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where K indicates the maximum connectivity. This value 
is shown in table I, where we can see that only totally 
unrealistic removals break the food web. 
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FIG. 4. Effects of random and directed removal on the con- 
nectivity distribution for the Silwood Park web for different 
fractions of removed species: / = 0.065 (10 removed species), 
and / = 0.1 (15 removed species). The estimated slopes are 
indicated ((a) r 2 = 0.87; p < 0.01; (b) r 2 = 0.87; p < 0.01; (c) 
r 2 = 0.92; p < 0.01; (d) r 2 = 0.77; p < 0.01. 



The nature of the behaviour of the three food webs 
is very similar despite differences in f c ,S and < s > 
between them. They exhibit high homeostasis when ran- 
dom species removals occurs, showing slow, linear de- 
crease in the fraction of species contained in the largest 
cluster S. The graph cannot be fragmented until ex- 
tremely high removal has been introduced. This can be 
seen in the values of < s >, which remain at (no species 
clusters different from S), or 1 (due to very few isolated 
species). But what is more revealing is that extinction 
rates remain at low values even for high /, so secondary 
extinctions are almost nonexistent. In fact we can esti- 
mate the fraction of removed species required in order to 
get food web fragmentation from random removal (Cohen 



However, what happens when most connected species 
are successively removed is clearly different. These webs 
are extremely vulnerable to such sort of removals. This 
fragility can be seen from: (a) the quick decay of S up 
to a critical fraction of removed species f c (see table I); 
(b) the high fragmentation of the food webs into species 
clusters disconnected among them, giving maximum lo- 
cal values of < s > at critical points f c where percolation 
takes place, and (c) the large fraction of isolated species 
occurring at low values of removed species, which reveal 
how fast secondary extinctions will occur (this is specially 
dramatic for our best defined web (Silwood), see figure 
3.e). A measure of this phenomenon is given by 7r c , the 
fraction of removed species that leads to an extinction 
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rate of one half. Although these are estimations based 
on a non-directed, non-weighted graph, other ingredients 
will presumably worsen this scenario due to other effects 
derived from indirect interactions (Yodzis, 1988; Stone 
and Roberts, 1991; Pimm, 1991) or habitat fragmenta- 
tion processes (Sole and Bascompte, 2001). 
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Table I: Summary of the average properties of the eco- 
logical networks analysed (Due to the irregular shape of 
P(k), no good power-law fit for Little Rock Lake can 
be obtained). (Here: (1) r 2 = 0.83; p < 0.01; (2) 
r 2 = 0.79;p < 0.01). 

We can make a simple division of the trophic nature of 
the species for each community into three groups: top 
predators, intermediate species and basal species (Po- 
lis, 1991). We find differences between analysed food 
webs in terms of the group that contains more highly- 
connected species, focusing in the set of species that are 
removed before reaching f c . For Ythan estuary, these 
are mainly intermediate species (fishes and invertebrate 
organisms, 60%), few top predators (birds, 20%) and 
parasites (15%) that cannot be easily included in any 
of these three groups. For Silwood park, most of these 
species are herbivores (66%), that could be considered 
as basal species since only one plant ( Cytisus scoparius) , 
is present. Hemipterous omnivores are also important 
(26%), but, as happens with parasites at Ythan, they do 
not belong to any of the three trophic groups. Finally, for 
Little Rock no basal nodes are highly connected, being 
intermediate species belonging to zooplankton, benthic 
invertebrates (70%) and top predators (such as fishes, 
24%) the most connected ones. 

Another fundamental question related to the fragility 
of these communities is how species removals affect the 
distribution of links per species under each kind of simu- 
lated attack. In other words, do food webs maintain their 
SF distributions as species are successively deleted? Our 
analysis shows that indeed, SF distributions are stable 
up to high values of / (more than 50% of species elim- 
inated in the three food webs) when random removals 
occur, while those particular distributions quickly disap- 
pear under directed removal. In Figure 4 we have repre- 
sented this for the Silwood web. Under random removal, 
the long-tailed distributions show little variation. In con- 
trast, that SF topology is lost when just few highly con- 
nected species are removed, which is likely to promote 
ecosystem collapse. 



IV. DISCUSSION 

The trophic organization of species-rich communities 
is similar to other complex network topologies (Albert et 
al 2000, Jeong et al 2000). They are extremely hetero- 
geneous, being their topology dominated by few highly 
connected nodes around which the rest of the network 
is organized, with a scale-free distribution of connec- 
tions. This complex organization entails some keys for 
ecological fragility. We have shown that SF food webs 
are very robust under random deletion of species. Sec- 
ondary extinctions remain at low values because the 
probability of being removed decreases with k accord- 
ing to P(k) w fc~ 7 , so it is unlikely that a highly con- 
nected species will be deleted. Such robustness becomes 
weakness under removals directed to species with many 
connections. We find specially significant the differences 
in extinction rates between random and directed attacks, 
raising up to 95 times higher for removals of highly con- 
nected (keystone) species. 

All the definitions of keystone species found in ecolog- 
ical literature share one feature: keystone species have 
large effects on other species in the community. The 
set of effects involved in each definition are very dif- 
ferent. Those effects have been studied mainly quali- 
tatively by the removal or introduction of species (sec 
Jordan ct al 1999 for a review). Quantitative approx- 
imations have used simulated food webs (Pimm 1991; 
Jordan et al 1999). In this respect, our approximation 
to the fragility of real, species-rich food webs through 
topological changes may help to design new quantitative 
methods for a priori identification of keystone species. 
We can identify keystone species as highly connected be- 
cause of the effects of their removal in terms of secondary 
extinctions. The quick arrival to f c and the high extinc- 
tion rates corresponding to low values of / when those 
keystone species are removed stresses the importance of 
identifying and protecting highly connected species that 
maintain the stability of ecological communities. 

By using this method, it is the topology of the food 
web instead of the trophic position of species what de- 
termines which species are keystone. In this respect, 
not only top predators must be considered as keystone 
species but also other organisms from different trophic 
levels, in agreement with previous studies (Bond 1993; 
Davie 2000). Making a simple division of species into the 
three classic trophic categories, we have seen that key- 
stone species belong to different categories in each of the 
analysed networks. A common feature found in the two 
best taxonomically described food webs is that species 
that feed on more than one trophic level (omnivorous 
species and parasites, which in most cases could be con- 
sidered as a special type of omnivory, Polis and Strong 
1996) are a representative group in the set of keystone 
species, enhancing the importance of omnivorous species 
in the stability of ecological communities (McCann and 
Hastings, 1997). 
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Our results suggest that there are basic principles of 
ecological organization (not revealed by previous analy- 
ses) underlying the assembly process of diverse communi- 
ties. These principles are in fact present at other scales: 
when a given spatial habitat is fragmented, there is a crit- 
ical percolation threshold which leads to the breakdown 
of the habitat available into many small patches (Bas- 
compte and Sole, 1996; Hanski, 1999). This threshold 
has very important consequences for metapopulation per- 
sistence and allows to define appropriate criteria for con- 
servation in fragmented landscapes (Keitt et al., 1997). 

Small-world and SF properties allows high food web 
complexity in terms of biodiversity. By one side, recent 
studies (Solow and Beet 1998; Montoya and Sole 2000) 
have shown that real food webs are more clustered from 
what would be expected from random wiring. These ev- 
idences support the hypothesis that compartmentaliza- 
tion (a characteristic of SW) is a way to enhance species 
coexistence in species-rich communities (May 1974). Be- 
sides, species that interact with a great number of species 
do so weekly. Conversely, species with strong interactions 
interact with few species of the community (May 1974, 
Paine 1992). Data on interaction strengths in natural 
food webs show that these networks are characterized by 
many weak interactions and very few strong interactions 
(Paine 1992; Raffaeli and Hall 1996). Recent studies also 
suggest that the balance of nature is related with those 
widespread weak interactions (McCann et al., 1998; Po- 
lis 1998). Thus highly connected species are important 
in promoting community stability and persistence. The 
detailed knowledge of this particular organization could 
lead to an improvement in the understanding of ecologi- 
cal assembly rules (Drake, 1990a; Drake 1990b). 

This approximation to community fragility has obvi- 
ous caveats derived from a lack of dynamics and taxo- 
nomic limitations. Besides, our findings are based on a 
non-directed, and non-weighted graph. Top-down and 
bottom-up effects due to species removals are considered 
together. A separate knowledge of these effects has been 
previously reported, arguing that whole community ef- 
fects were more interesting and more relevant in the de- 
termination of keystone species (Jordan et al 1999). Pre- 
liminary analysis of model ecosystems suggest that our 
results are robust (Sole and Montoya, unpublished). 

Current estimations of extinction rates are based on 
species-area relations, combined with estimates of habi- 
tat loss (May et al., 1995). The addition of secondary 
extinctions due to removal of keystone species, together 
with other indirect effects are likely to increase such pro- 
jections. Food webs are described at a local scale, but the 
estimated extinction rates obtained from our study could 
be important in forecasting extinction rates at regional 
and global scales. 
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